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Determining the equation of state of dark energy with astronomical observations is crucially im- 
portant to understand the nature of dark energy. In performing a likelihood analysis of the data, 
especially of the cosmic microwave background and large scale structure data the dark energy per- 
turbations have to be taken into account both for theoretical consistency and for numerical accuracy. 
Usually, one assumes in the global fitting analysis that the dark energy perturbations are adiabatic. 
In this paper, we study the dark energy isocurvature perturbation analytically and discuss its impli- 
cations for the cosmic microwave background radiation and large scale structure. Furthermore, with 
the current astronomical observational data and by employing Markov Chain Monte Carlo method, 
we perform a global analysis of cosmological parameters assuming general initial conditions for the 
dark energy perturbations. The results show that the dark energy isocurvature perturbations are 
very weakly constrained and that purely adiabatic initial conditions are consistent with the data. 



I. INTRODUCTION 



Since the discovery of the accelerated expansion of the universe by observations of distant Type la supernovae 
(SNIa) in 1998 (T), [||, dark energy has become a hot topic in physics and astronomy. So far a lot of models have 
been proposed in the literature. In general these models can be classified according to the equation of state (EoS) 
w e of the dark energy, defined as the ratio of its pressure to energy density. The simplest assumption is to consider 
dark energy with constant w e , and more specifically to assume a cosmological constant whose EoS is equals to — 1. 
Although this scenario is consistent with observations it suffers from the well-known fine-tuning and coincidence 



problems! 



|. Alternatively, dynamical dark energy models, such as quintessence phantom k-essence fl"o| . 

quintom [111 |12| and so on, have a time-dependent EoS. For quintessence w e > —1, while for phantom w e < —1. But 
for quintom models, the EoS crosses the boundary set by w e = — 1. To investigate dark energy without making use 
of specific field models, one often parameterizes the EoS of dark energy asfl3j 

w e (a) = w + w a (l - a) (1) 

where a is the scale factor normalized to be 1 at the present time. We adopt this parametrization in the current work. 

Given the fact that a lot of theoretical models exist in the literature, it is crucially important to use the accumulated 
high precision observational data from SNIa, the cosmic microwave background (CMB) radiation and large scale 
structure (LSS) surveys to constrain the value of w e (a). Since dynamical dark energy should fluctuate in space as 
described by the conservation of the energy-momentum tensor, in a global analysis with a general time evolving 
EoS one should take into account the dark energy perturbations in order to have a consistent procedure. This is 
particularly important when fitting cosmological parameters to the data of CMB and LSS. Simply switching off the 
dark energy perturbation is not theoretically correct and will lead to biased results. Numerically it has been shown 
that the results obtained are quite different between the two cases with and without the dark energy perturbations 

[Mill. 

With the generally parameterized EoS, there inevitably exists a singularity when w e = — 1. When w e crosses this 
critical point, the dark energy perturbations will diverge [111 [22H24| . It has been shown that in the context of general 
relativity, it is impossible to obtain a background which crosses the "cosmological constant boundary" with only a 
single scalar field or a single perfect fluid. In fact, this is the reason why the quintom scenario of dark energy needs 
to introduce extra degrees of freedom [ll|, US 122-28] l . In order to handle the perturbation when w e crosses — 1, 
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1 For a consistent and complete proof of the no-go theorem, please see, Il9ll . 
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we proposed a method to deal with the dark energy perturbations during the crossing of the boundary w e = — 1 in 
Ref. [20l | . According to this method, the energy and momentum density perturbations of dark energy are treated as 
constant during the small interval around the critical point w e = — 1. This method is justified in Ref. [2!| from the 
viewpoint of general relativistic matching conditions. 

Another issue concerned with the dark energy perturbations is the question of initial condition. Generally, there are 
two types of initial conditions for the perturbations: adiabatic and isocurvature. In Refs. [H], Ell for the global 

fitting of the dark energy EoS to the observational data, it was assumed that the perturbations were purely adiabatic. 
In this paper we will study more general initial conditions which admit dark energy isocurvature perturbations and 
discuss the implications for CMB temperature and polarization power spectra and the LSS matter power spectrum. In 
the literature, baryon and dark matter isocurvature perturbations have been extensively discussed and tight constraints 
on these are obtained. There have also been studies of the dark energy isocurvature perturbations which, however, 
are usually limited in the framework of quintessence models [30l - [34l ]. In these studies, it has been shown that the 
quintessence isocurvature perturbations could lead to the suppression of the CMB quadrupole via the anti-correlation 
between the adiabatic and the isocurvature modes (3314361 ]. In this paper, working with the parameterized EoS, we 
consider both adiabatic and isocurvature initial conditions and the correlation between them in the likelihood data 
fitting analysis. We discuss the current constraints on the cosmological parameters, with result when admitting the 
possible existence of dark energy isocurvature modes. Our paper is organized as follows: in section II, we briefly review 
the theory of perturbations; In section III, we analytically study in detail the dark energy isocurvature perturbations; 
In section IV we study effects of the dark energy isocurvature perturbations on CMB and LSS, and we present the 
current constraints on them in Section V; Section VI is our summary. 



II. ADIABATIC AND ISOCURVATURE PERTURBATIONS 



We consider a spatially flat Fricdmann-Robertson- Walker universe as the background. The metric of the perturbed 
spacetime in the conformal Newtonian gauge reads, 

ds 2 = a{T]) 2 [(l + 2$)drj 2 - (1 - 2$)5 ij dx i dx j ] , (2) 

where we have implicitly assumed that the shear perturbations can be neglected and the metric perturbations are 
fully described by one relativistic potential $. In the matter sector, the perturbations arc expressed by the perturbed 
energy-momentum tensor which is gauge dependent. However, for the discussions of perturbations on large scales it 
is more convenient to use gauge invariant variables constructed by combining the energy-momentum perturbations 
with the metric perturbations. In this paper, we use the following gauge-independent variables for each species 

U 3(1 + ' 

A Q = + -j^(Pa +Pa)0 a , (3) 

where S a = Sp a /p a is the density contrast, 6 a = ik l ST° ia /(p a + p a ) is the corresponding momentum density per- 
turbation, and the conformal Hubble parameter is defined by % = a' /a with the prime denoting the derivative with 
respect to conformal time. £ a is a comoving curvature perturbation, and as we will see later in this paper A Q may be 
called an "effective" density perturbation. The conservation of the energy-momentum tensor at the linear order gives 
the equations governing the evolutions of £ a and A a : 

C + 3W(cJ a - <4)-^L_ + ^(-^- Ca) = |U , (4) 

Pa + Pa 37t Pa + Pa 3% 

K + + ^)A a - (U - ?L + ^-)(p a+PaKa = (p a +p a )[*> + (2H ^ + ^)$] . (5) 

In the above equations, c sa is the sound speed defined in the comoving frame of the fluid while the so-called adiabatic 
sound speed, c aa , is defined as c 2 a = p' a / p' a = w a — w' a /[?>H{l + w a )}. For a perfect fluid c sa = c aa , and for a 
canonical scalar field c sa = 1. 

To close the system, we also need the Poisson equation, 

^0> = -12^?^A Q , (6) 
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which can be obtained from the perturbed Einstein equations. Comparing this equation with the Poisson equation in 
the Newtonian gravity we can see that A Q may be called an effective density perturbation. On super horizon scales 
k?] <C 1, the terms proportional to fc 2 in Eqs. flU, © and (jH) can be dropped, and these equations become 

C+3H(c 2 sa -cl a )-^-=0, (7) 

Pa + Pa 

K + m - ^)A Q -{%- ^)( Pa+Pa )C a = {p a+Pa )W + (2H - ^)$] , (8) 
0. (9) 



From equation ([7]), we see that for a perfect fluid £ a is conserved on large scales. 

To solve the set of perturbation equations, we need to specify the initial conditions. Usually these initial conditions 
are set at the time deep inside the radiation dominated era when the scales corresponding to observations today were 
far outside the horizon. Inflation provides a natural mechanism to generate these initial perturbations. They originate 
from quantum vacuum fluctuations and become classical perturbations when their corresponding length scales leave 
the horizon during inflation. In the post-inflation epoch, these primordial perturbations re-enter the horizon and 
interact with matter to cause the CMB anisotropies and structure formation. Hence it is important to inspect the 
evolutions of the perturbations on super-horizon scales until the scales re-enter the horizon. For this purpose we only 
need to consider the equations (J7J, © and ©• There are two types of solutions of these equations, called adiabatic 
and isocurvature (or entropy) modes. For adiabatic perturbation all comoving curvature perturbations £ a are the 
same as that of radiation £ r . The total comoving curvature perturbation is also equal to £ r , 

C= — ~ V(/5 Q +p«)Ca =Cr , (10) 
n -\- m — ' 



which is constant since C' r = 0- Thus, for adiabatic perturbation, the picture is simple: the primordial perturbations 
C are frozen while they are outside the horizon. With Eq. ©, the sum of Eq. © over all species gives 



Integration of this equation gives 



(^-•H)C = *' + (2«-^)*. (11) 



^ = c£ f^r), (12) 

a z a A J H(a) 



where C is a constant. The first term on the right hand side decays in the expanding universe and can be neglected. 

However, if one of the comoving curvature perturbations £ a is not equal to £ r , the density perturbation have an 
isocurvature mode, the isocurvature perturbation of species a is defined by 

S a = 3(Ca ~ Cr) = ~t~T~ - l S r ■ (13) 

1 + w a 4 

In principle, if the universe contains TV components, there should be at most N — 1 isocurvature density perturbations. 
With the presence of isocurvature modes, the equation (jlip relating the potential and the total comoving curvature 
perturbation is still valid. However in this case C 7^ Cr- If only one species a has isocurvature perturbation, £ is 

^ + 5&T3 S «' <»» 

It is not conserved on large scales. We may define £ Q = (p a + p a )S a /3, so that 

Ja_ 

P + P 



C = (r + —, (15) 



and Eqs. ([7]) and (J8j) become 

C + 3H{1 + cL)£a + 3H(c 2 Q - cL)A« = 



k + m ^)a q + &-u)(i- p -^^)a a - . (i6) 

n n P + P 
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The above equations describe how the isocurvature perturbations evolve on large scales and Eq. (|T4|) or (fT5|) char- 
acterizes the contribution of isocurvature perturbations to the total comoving curvature perturbation. The potential 
can be solved by integrating Eq. (JTTJ) and we obtain 




where $ adl is the contribution of the adiabatic mode given in Eq. (|12[) . and the last term is the contribution from 
the isocurvature perturbation. To get it we have used the equation W — Ti 2 = —Ai:Ga 2 {p + p). This equation shows 
explicitly that both adiabatic and isocurvature perturbations are able to generate the metric perturbation. 

If the species a is subdominant (p a + p a )/(p +f>) 0, its perturbation has a negligible contribution to the metric 
perturbation and the potential $ ~ <J> adl . Then the evolution equations (fTH]) of isocurvature perturbations become 

C + m{\ + cL)£ Q + 3H(cL - cL)A« = o , 

K + m-^)A a + (^-H)Z a =Q . (18) 

These two equations are the basis for the discussion of dark energy isocurvature perturbations in the radiation and 
matter dominated eras in the next section. We can see from Eq. (|15j) that the contribution of dark energy isocurvature 
perturbations relies on the ratio Ce/(p+p) compared with £ r , where the subscript e represents dark energy. Because 
Cr is conserved, qualitatively the effect of the isocurvature depends on whether £ e /(p + p) grows or decays with 
time. When the density of dark energy becomes significant at late time, its isocurvature perturbations could make as 
important contribution to the metric perturbation and we should use the equations (fTH]) to investigate its evolution. 



III. DARK ENERGY ISOCURVATURE PERTURBATION 



In this section we will study the dark energy isocurvature perturbations during the radiation and matter dominated 
epochs. For simplicity we assume that the perturbations of baryons, dark matter, neutrinos and so on are adiabatic. 
The dark energy was subdominant in the early universe and it only starts to dominate the universe at low redshifts. 
So we can use the equations (|18jl to study the dark energy isocurvature perturbation on super-horizon scales. 

First of all, we will discuss the behavior of dark energy isocurvature perturbations for some specific dark energy 
models. 



A. Single Fluid 

If dark energy is a perfect fluid with one component, its sound speed in its comoving frame is equal to the adiabatic 
sound speed, c 2 se = c 2 e . It cannot be negative otherwise its perturbation would be unstable on small scales. From the 
first equation of ([Tg)) . we have 

£ + 3«(1 + = . (19) 

This equation implies 

£ e OC p e +p e . (20) 

So, the contribution £ e /(p + p) of dark energy isocurvature perturbation scales as (p e + p e )/(p + p). If c 2 e changes 
slowly, we can treat it as a constant. In this case, £ e /(p + p) will not be damped if c 2 e < 1/3 in the radiation epoch. 
However, £ e /(p + p) will decay in the matter dominated era unless c^ e = 0. 

B. Multiple fluids 

If the dark energy contains multiple perfect fluids, then for each component we have £ e i = Ci(p e i +p e i) and 

£,e = = ^2 C i(Pei + Pei) , 

i i 

= T,i C i(Pei + Pei) ^ 
P+P P+P 
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where Ci are constants. If all the constants Ci are the same, then there are no internal isocurvature perturbations 
among the components of dark energy, £ e /(p+f>) scales as (p e + p e )/(p + p) like in the case of single fluid. Otherwise, 
we should study the evolution of the isocurvature perturbation for each component individually. 



C. Single Field 



There are many proposals for dark energy models based on scalar fields. Simple dark energy models like quintessence, 
phantom and k-essence can be constructed from one scalar field. For quintessence or phantom dark energy, the sound 
speed is c 2 e = 1. For k-essence, c 2 e can have any nonnegative value. In general the behavior of scalar field dark energy 
models cannot be solved analytically. But in some cases, c 2 e and the equation of state w e change very slowly and can 
be treated as constants, and consequently the analysis becomes simpler. In this case the equations (fl~8|) reduce to the 
following second order differential equation for £ e , 

C + [(7 + 3w e )H - 2%-]£ + 3[(4 + 3w e + c 2 se )H 2 - (1 + c 2 e )U'% = . (22) 
H 

In the radiation dominated era, a oc 77 and "H = I/77, and the solution of Eq. (|22ll is 

£ £ = r? -(8+3»e)/2 [Cir? |V( 3w =+ 2 ) 2 - 24c ia + C 2V -^ 3W ° +2 ^- 24C "} . (23) 

The contributions of dark energy isocurvature perturbations to the metric perturbation scales as 

— ^_ OC ,,-3™./»[Ci^V(3».+ 2 > a -Mc». + -VC*».+2) a -2**.] . (24) 
Pr + Pr 

For quintessence or phantom dark energy, c 2 e = 1 and w e < 1, we can see from the above equation that the 
effect of the isocurvature perturbations of any model with positive w e decays with the expansion of the universe. If 
w e = 0, £, e /(Pr + Pr) oscillates with constant amplitude. For negative equation of state, the effect of the isocurvature 
perturbations grows with time. An interesting case is that when the scalar field is almost frozen during the radiation 
epoch, i.e., w e ~ —1, £ e /(p r + Pr) grows as a 3 / 2 even though £ e and S e decrease with time. 

For k-cssence dark energy, c 2 e may not be equal to one, but it should be nonnegative. In this case, for the models 
with smaller sound speed, the contribution of its isocurvature perturbation is more likely to be growing. 

Similarly in the matter dominated era, a oc rj 2 and TL = 2/77. Thus, for constant w e and c 2 e , one has 



P+P Pra 

This shows that the contribution of dark energy isocurvature perturbations has a similar behavior as in the case 
of radiation dominated era as discussed above. With more negative equation of state or smaller sound speed, this 
contribution is more likely to grow with time. 



D. Multiple fields 

If dark energy contains multi-fields, like in the case of the quintom model, we should solve the equations (|18[) for 
each component. Here, for simplicity, we have assumed that there are no interactions other than gravity among the 
internal components of the dark energy. The total contribution of dark energy isocurvature perturbations is 

6 _ Ei U (26) 



P+P P+P 

For example, consider the quintom model with two fields, one quintessence and the other a phantom field, and 
assume that each field has an extremely small mass. Both fields are slowly rolling in the radiation dominated era and 
£,ei/(pr + Pr) oc a 3 / 2 for i = l, 2. Hence we have 

ie oc a 3 / 2 . (27) 



Pr + Pr 
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IV. THE EFFECTS OF DARK ENERGY ISOCURVATURE PERTURBATION 



lJ ~ \ VA adi A S oCosA A iso I ' 



In this section we study the effects of the dark energy isocurvature perturbations on CMB and LSS observations. We 
will take the paramctrization w e (a) = w + w a (l — a) and consider the sound speed c 2 e as an arbitrarily non- negative 
parameter. Besides, we also need to parameterize the power spectra of the initial perturbations. Statistically, both 
adiabatic perturbation £ r and isocurvature perturbation S e are treated as random fields as predicted by inflation 
theory. To be general, we should consider the correlation between them. If their statistics are Gaussian, both the 
adiabatic and isocurvature fields arc fully described by the power spectra. To characterize a well defined system 
including both adiabatic and isocurvature modes, one usually introduce a vector Afj with two components, 

x = jCr adiabatic , 

I S e isocurvature . 

Then the primordial power spectra Vij are defined by 

2tt 2 

(Ai(k)Ar;(k')} = -p-P«(*)«(k - k'). (29) 

One can parameterize the power spectra as Vij = A i j(-£-) ni '~ 1 , where and are 2— dimensional matrices 
which characterize the amplitudes and spectral indices, respectively. We have 

^adi \MadiAso COS A 

where cos A = -^=s|= describes the correlation between adiabatic and isocurvature perturbations [37}, and A ad i 
and A iso are the amplitudes of adiabatic and isocurvature modes respectively. The spectral indices are denoted by 
n^ dl = nn and nf° — n 2 2- For simplicity we assume that n c s m ' = n i2 = Ttll +" 22 [38| . 

As we see from Eq. (|17|) . both adiabatic and isocurvature modes can generate metric perturbations and therefore 
temperature anisotropics. Symbolically we have 

ST \ Cu\ 

~7^~ — ("^"Jadi + ("^Tjiso, (31) 

and hence, the temperature angular spectrum C% can be expressed as 

Ci = A adi C? di + A iso Cl so + 2v/A adi A iso cos ACf di ' iso (32) 

where 

4-7T f k 

C? = J d\nk(-r^Ql(k)el(k) (33) 

with ®\ being the transfer function of photons for the initial condition i. There are similar formulas for the CMB EE 
and BB polarization spectra and temperature-polarization spectrum TE. 

The isocurvature perturbations also affects the matter power spectrum P{k) as follows, 

P(k) = A adi P* di (k) + A iso P iso (k) + 2^/ A adi A iso cos AP adi > iso (fc) , (34) 

where P 1 * (k) can be described as 

ptf(fc) = (^) n *- x T i (k)THk), (35) 
k 

with T' l (k) being the transfer functions of matter perturbation for initial condition i. 

In order to show the effects of the isocurvature perturbations on CMB and LSS observations, we plot in Fig|T]the TT 
and TE power spectra of CMB and in Fig |2] the matter power spectrum in the case of fully anti-correlation, cos A = — 1. 
In the computations, the fiducial cosmological parameters are chosen as wo = —1.148, w a = 1.01, c 2 e = 0.01, 
A adi = 2.36 x 10" 9 , nf l = 0.95, A iso = 1.48 x 10~ 8 , nf° = -1.7, tu b = 0.02247, uj c = 0.1135, H = 71.8km/s/Mpc, 
where ojj, = Q,f,h 2 and ui c = Vl^ti 2 denote the physical baryon and cold dark matter density parameters, respectively, 
and Hq = 100/ikm/s/Mpc is the current Hubble constant. The main effects of isocurvature perturbations of dark 
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FIG. 1: Top Panel: The angular power spectrum of CMB. Bottom Panel: The TE power spectrum of CMB. The red solid lines 
denote the spectrum obtained including the contribution of anti-correlated adiabatic and isocurvature perturbation, while the 
black dashed line is obtained by only including the adiabatic contribution. 
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FIG. 2: The matter power spectrum obtained with the cosmological parameters chosen to be the same as in Fig[T] The red 
line denotes the total power spectrum and the black dash-dotted line is that with only the adiabatic component. 

energy are on large scales. One can see the suppression in the CMB quadrupole which is realized by the anti-correlation 
between isocurvature and adiabatic perturbations. However, in the matter power spectrum, there is an increment on 
large scales with this set of parameters chosen. 

Since the effects of isocurvature perturbations appear mainly at large scale, it will be difficult to get a tight constraint 
on it with current data. This is because we know that in the case of the CMB, the data on large scale is cosmic 
variance uncertainty dominated, while for LSS the largest scale we can observed today is only about k = 0.02/iMpc _1 
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V. OBSERVATIONAL CONSTRAINTS 
A. Data and Cosmological Parameters 

We extended the publicly available MCMC package CosmoMC 2 by including dynamical dark energy and its 
perturbations discussed in this paper. We then performed a global analysis. In the computation of the CMB we 
have included the WMAP7 temperature and polarization power spectra with the routine for computing the likelihood 
supplied b y th e WMAP team 3 141| . Furthermore, we include small scale temperature anisotropies measured by ACB AR 
[42 ] , CBI [43j and Boomerang [441 ] . The matter power spectrum measured by observations of luminous red galaxies 
(LRG) from SDSS [3^], and the "Union II" supernovae dataset [4^| was also taken into account. Furthermore, we 
added a prior on the Hubble constant, Hq = 74.2 ± 3.6 km/s/Mpc given by ref. 14611 as well as a weak Gaussian prior 
on the baryon density ojb = 0.022 ± 0.002(lcr) from Big Bang Nucleosynthesis [47j 4 . Simultaneously we also used a 
cosmic age tophat prior as 10 Gyr < to < 20 Gyr. 

In the numerical calculation we considered the most general parameter space 

P = w c , G s , t, w , w„, c^ e , cos A, nf , 4 di , 4 } , (36) 

where S = 100^- is the ratio of the sound horizon to the angular diameter distance at decoupling and r characterizes 
the optical depth to reionization. 

B. Global Fitting Results 

In order to show explicitly the effect of dark energy isocurvature perturbation, we have done two different kinds 
of calculations: one is with pure adiabatic initial condition by making the three parameters in Eq. ([36| . i.e. nf°, A; so 
and cos A vanish and another is the calculation with the full set of parameters in Eq. (|36p including the dark energy 
isocurvature perturbation (hereafter refer to as "Mixed"). First of all we considered the effect of the isocurvature 
perturbation on the determination of the dark energy EoS. In TabU we present the numerical values of the dark 
energy EoS. 





too 


1t)a 


Adiabatic 


-1.143 ±0.160 


0.463 ± 0.605 


Mixed 


1 ,nr, + 0.1B4 

i.ioz_ 158 


n41 n + U.BU2 



TABLE I: Current limits on the dark energy EoS 

One can see from this table that, given the current observation the isocurvature perturbations makes a change. 
However the effect is small. As one can see from FigfT] the cause of the change is the isocurvature mode which 
suppresses the power spectrum on large scale. Theoretically, to have a sizable suppression, a smaller w e (a) is required 
in the early universe. This explains why a smaller mean value of w e (a) is obtained in the mixed case given the current 
data with the suppressed CMB quadrupole. In Fig [3] we plot the two dimensional constraints on the dark energy EoS 
parameters wq, w a . To show the importance of dark energy perturbations we also present the results obtained when 
the DE perturbations switched off incorrectly We can see from this plot that it brings an error 9% on wo and 52.1% 
on Wa- 
in FigUl we plot the marginalized 1 — D probability distribution of the parameters related to the initial conditions. 
The constraints on isocurvature parameters, such as Ai SO , nf° and cos A of dark energy are weak. This is understand- 
able, since, as shown in FigfT] and [2J the isocurvature perturbations of dark energy mainly make contribution at large 
scale, where the observational data are limited. We note that for the mixed case, the mean value of the adiabatic 
primordial perturbation amplitude A ac ji ~ 2.5 x 10 -9 is slightly larger than in the adiabatic case A a ^ ~ 2.4 x 10 -9 , 
while the spectral index is smaller. This indicates that, with respect to adiabatic perturbation, the isocurvature mode 
has a negative effect on the power spectra on large scale (small k), i.e., the mixed mode can suppress the CMB TT 
angular power spectra at low values of I. Moreover, the decrease of x 2 , A% 2 = xidi ~ Xmix ~ 2 gives a hint that the 
mixed case is mildly favored by the data. 



2 http : //cosmologist . inf o/cosmomc/ . 

3 Available at the LAMBDA website: http://lambda.gsfc.nasa.gov/. 

4 Actually, the chosen of prior may lead some uncertainty. However, we have checked that the change is neglectable. 
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FIG. 3: The constraints on the dark energy EoS parameter Wo and w a . The red solid lines are for the mixed case and black 
dash-dotted line for adiabatic case. The blue dashed lines stand for the constraints without dark energy perturbation. 
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FIG. 4: The marginalized ID probability distribution of the cosmological parameters related to dark energy perturbation. 
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VI. SUMMARY AND DISCUSSION 

As is well known, cosmological perturbations are crucially important for understanding the CMB anisotropics and 
structure formation. Up to now the theory of cosmological perturbation is very successful and has been confirmed 
by the high precision observations and experiments such as WMAP and SDSS. Since the perturbed spacetime is 
determined by the perturbations of all of the matter components in the universe, it is also important to study the 
dark energy perturbation. If we naively switch off the dark energy perturbation, the result would be misleading. 
Moreover, to be general, besides the adiabatic perturbation which is mostly studied in the literature, one should also 
consider the isocurvature perturbation. Because dark energy couples very weakly to other matter it is not so easy 
to construct a dark energy model which has purely adiabatic perturbation. Dark energy isocurvature perturbations 
have important application to lower the quadrupole of the CMB angular power spectrum as needed by the COBE 
and WMAP observations^. 

In this paper, we have studied in detail the effects of dark energy isocurvature perturbations. We have included 
dark energy isocurvature perturbation in the data analysis. By employing a Markov Chain Monte Carlo method, 
we have performed a global analysis of the determination of the cosmological parameters from current astronomical 
observational data. We find that isocurvature perturbations decease x 2 by about 2, and has small effect on other 
parameters. The current limit on the isocurvature initial condition is weak. We expect that future precision measure- 
ments of CMB and LSS on large angular scales, especially the measurements of CMB-LSS cross correlations will lead 
to a tighter constraint. 
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